A gauge-invariant and current-continuous microscopic ac quantum transport theory 
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There had been consensus on what the accurate ac quantum transport theory was until some 
recent works challenged the conventional wisdom. Basing on the non-equilibrium Green's function 
formalism for time-dependent quantum transport, we derive an expression for the dynamic admit- 
tance that satisfies gauge invariance and current continuity, and clarify the key concept in the field. 
The validity of our now formalism is verified by first-principles calculation of the transient current 
through a carbon-nanotube-based device under the time-dependent bias voltage. Moreover, the 
previously well-accepted expression for dynamic admittance is recovered only when the device is a 
perfect conductor at a specific potential. 
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Understanding quantum transport is important for the 
development of nanoelectronics. The Laudaure-Biittiker 
formula has been widely used to calculate steady state 
current with great success. [THS] The resulting steady 
state currents through the leads in contact with the elec- 
tronic device satisfy gauge invariance and current conser- 
vation. Here, gauge invariance condition means that all 
physical observables should remain unchanged when the 
applied bias voltages are shifted by a constant. Current 
conservation implies that the steady state currents sat- 
isfy the Kirchhoff's circuit law, i.e., the currents through 
all leads sum up to zero. In contrast, the situation is 
somewhat different for ac quantum transport. The gauge 
invariance should still hold. Meanwhile, the particle cur- 
rents through all leads are no longer conserved, due to 
the possible transient charge accumulation or depletion 
at the device. In such cases, the current conservation 
(or Kirchhoff's law) does not apply. Instead, the cur- 
rents satisfy the charge continuity equation. Of course, 
the current conservation may be retrieved formally by 
introducing the displacement current, so that the total 
(particle plus displacement) currents through all leads 
sum up to zero. 

In 1993, Biittiker and coworkers have derived an ex- 
pression for the linear dynamic admittance of a small 
conductor by employing the scattering matrix theory. [7] 
By taking into account of a uniform potential distribution 
inside the small conductor, gauge invariance was guar- 
anteed. To ensure current conservation, the conductor 
was treated as an extra terminal, and the resulting ad- 
ditional dynamic admittance matrix counts for the con- 
tribution from the displacement current. However, it is 
unclear how to apply such a formulation to a generic de- 
vice where the potential distribution is not uniform. In 
1995, Anantram and Datta pioneered a microscopic ac 
transport theory based on nonequilibrium Green's func- 
tion (NEGF) formalism and the wide-band limit (WBL) 
approximation. f8| However, the resulting formalism for 
the dynamic admittance does not satisfy gauge invari- 
ance and current conservation. To remedy such a situ- 
ation, in 1999, a phenomenological NEGF method was 



proposed to ensure both gauge invariance and current 
conservation by Wang, Wang and Guo.[5j The key idea 
was to partition the displacement current onto each elec- 
trode or terminal in a phenomenological manner. This 
was achieved by imposing the requirements of both gauge 
invariance and current conservation on top of the theory 
of Anantram and Datta. |8j It was a very nice idea and 
the resulting expression for the dynamic admittance has 
been widely used in the field ever since, [TU| despite of the 
two facts: (1) it was built on a previous formalism that 
is not gauge-invariant; [8] and (2) no detailed justification 
or derivation based on the microscopic theory was given. 

We calculated the dynamic admittance of a carbon 
nanotube (CNT) based device, and compared it with the 
result of a well-established NEGF-based time-dependent 
density-functional theory (TDDFT-NEGF) method. dH 
IT2| The two sets of results are found significantly dif- 
ferent from each other (see Fig. lb). Recently, Wang 
and coworkers have calculated the dynamic admittances 
of a benzene-dithiol and a chain of carbon atoms sand- 
wiched between two aluminum electrodes. fl3j Both the 
phenomenological formalism of Reference [S] and a micro- 
scopic theory[T3] were employed in their calculations. It 
was found that the discrepancy between the two could be 
as much as two orders of magnitude. The inconsistency 
in the calculation results thus presents a serious challenge 
for the existing ac quantum transport theory. 

Before we proceed to develop a correct microscopic the- 
ory for ac quantum transport, we explore what it might 
be. Current conservation for steady state current comes 
from the current continuity equation V ■ J + — 0, 
where J(r, t) is the current density function and p(r, t) 
the electron density function. For steady states, the time- 
partial derivative of p{r,t) is zero, so V • J = 0, which 
leads to current conservation. However, for ac transport, 
p{f, t) is time-dependent and its time-derivative is usually 
nonzero, and thus V • J 7^ 0, i.e. the current conservation 
may not hold. In other words, for ac or time-dependent 
current, the current conservation is not required. Cer- 
tainly, if the displacement current is included, current 
conservation is satisfied for the total current, just as that 
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in the conventional theory of electrodynamics. 

As demonstrated by Biittiker and coworkers, gauge in- 
variance can be preserved by taking into account of the 
induced potential distribution of the electronic device 
properly.^ Due to the bias voltages on the electrodes, 
the electrons in the device region respond via redistri- 
bution. Moreover, electrons may leave from or enter 
into the device, which leads to further charge redistri- 
bution. The Hartree potential in the device then ad- 
justs to the applied voltages and charge redistribution. 
As a result, the relative Hartree potential distribution 
is gauge-invariant, so is the resulting current, as required 
by physics. Therefore, the key to ensure gauge invariance 
is to self-consistently determine the Hartree potential. 

Our objective is thus to develop an ac quantum 



transport theory that is gauge-invariant and current- 
continuous (rather than current-conserved). Starting 
from the Keldysh NEGF formula in time domain, [TS] 

IJt)^-J^ dtirr[G<(t,ti)S^(ii,i) 
+ G'^(i,ti)S<(ti,t) + /i.c.], 

we perform a double-time Fourier transform to convert 
it into frequency domain. The Green's function and self- 
energies are separated into equilibrium and small-signal 
terms as A'^ = + SA'', where ^ = r,a, <, or >, A = 
G or E. The frequency-dependent current can thus be 
expressed as 



Uco) ^'-j^ —Tr[G^,{E+)SJ:<{E+,E) + G<{E+)SJ:^^{E+,E) + 5G<{E+,E)j:UE) 

(2) 



where E^ = E + hoj. Employing Langreth's rules, [S] we 
obtain 



6G<{E,,E2) ^ Gl,{E,)S^<iE,,E2)G^oiE2) 
+ 6G^iE,,E2)^^iE2)G^{E2) 
+ G^oiEi)^o{Ei)SG''iE,,E2). 



(3) 



And the small-signal terms of the retarded and lesser self- 
energy are obtained by performing Fourier transform on 
time-domain expressions! 16) 
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where Uq. is the applied voltage at Lead a, Ta 
is the coupling between device and Lead a, and 
equilibrium self-energies are T,q^{E) — ~iTa{E)/2, 
E<^{E) = if{E)Ta,{E). The retarded Green's func- 
tion is expressed by the Dyson equation, G^{t,t') = 

g''{t,t')+ /dtidt2G'^(t,ii)E'-(ti,f2)G'-(t2,i'), where g"- 
is the Green's function for uncoupled the system, and 
its corresponding small-signal term Sg^ is the change of 
g^ due to the change of Hartree potential. Although un- 
coupled, the device still experiences a potential shift after 
the application of the external potential; therefore, dg^ is 
nonzero in general, which was neglected in Reference [8]. 



Therefore, 

SG^iE,,E2) - G^aiEi)MiEi)r'Sg^{E,.E2){g^oiE2)r' 
+ 5^-'\E^,E2)]Gl{E2), 

(5) 



and (5G" is the Hermitian conjugate of JG*" 
evaluated from 



can be 



Sg'-itiM) 



dcj 



U{uj){e 



'-)g^oiti-t2), (6) 



where U{f,uj) is the induced Hartree potential energy 
in the device region due to the bias voltages, which can 
be obtained self-consistently by solving Poisson equation 
subject to the proper boundary condition. fl7J Within the 
linear-response regime, U{r,uj) can be approximated as 
U{r,u!) = — eu^(r, w)ti.^(w), where Uj{f,t) is the di- 
mensionless potential distribution with boundary condi- 
1, and = 0, (7 9^ T,^Uj{f,uj) = 1 



tion V. 
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is required by gauge invariance. By performing Fourier 
transform, we have 



6g^{E+,E)=J2 



2hu} 



■K(i?)-55(£;+)]. (7) 



As a consequence, 6G^ can be expressed as 

SG^iE\E)^Y.^-I^G^oiE^) 
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r^(i?)-r^(i?+) 



(8) 



Substituting Eqs. ([3|, Q and (|8| into Eq. we obtain 
the expression of dynamic admittance for the particle 
current as: 



3 



_ ie^ r dE 

+ /G5 {u^(f,c^);ic^ ~ ^A_^} GS(r„ + ^^GgA+„) - /Gg {up{T, uj)huj + ik^p} GS(f „ - zA+^G^f )], 
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where for simplicity, A = A{E),A = A{E+),K^.y = 

T-,(E)^T-,(E+) . _ T^(E)+T^{E+) 
1 ^-1-^ — 2 



2 '^'-+7 



r = EaTa, and 
Eq. (9| is the central equa- 
tion of our microscopic theory. The gauge is naturally 
satisfied as Ga-y (oj) = based on Eq. (joj) . In gen- 
eral, current conservation is not satisfied for ac currents 
through a nano-device; instead, current continuity holds, 
which can be verified as J^-^ ^^^{lo) = i^Q{i^)\v^=id,i3^K., 
and 



Q(w) 



2uj 



+00 



dE 
2^ 



Tr[5G<{E+,E)] (10) 



is the accumulated charge of the device, and is gener- 
ally nonzero. We have thus developed a gauge-invariant 
and current-continuous NEGF formalism for ac quantum 
transport. In practice, we need to solve u^(f, w) self- 
consistently subjected to the boundary condition = 1, 
and Vfj = 0|^^^. When the frequency uj is low, the WBL 
approximation yields accurate dynamic admittance. Un- 
der the WBL approximation, Eq. ^ can be greatly sim- 
plified as 



+ 00 



AE 

hid 



rr[(r„GSrGS - iSt^r„GSGS)5, 



a/3 



TaGSr/^Gg + ihu;Glu^{r,u)G^oT^]if - J). 



(11) 



Eq. (11) is consistent with the formula for the dynamic 
admittance of the total current in Reference [T3]. The 
last term in Eq. ( 11 ) that is proportional to up{r, u) was 
mistakenly termed as displacement current, and the rest 
as the particle current. In fact, all terms are for the 
particle current. 

According to Reference [9|, within the WBL approxi- 
mation, the dynamic admittance can be expressed as 
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rr[(r„GsrGS - i;ic^r„G5GS)<5, 
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F^GSF^Gg + zA„;iwG5GgF^](/ - /), 



(12) 



where the first three terms are the results of Anantram 
and Datta,[8 and the fourth is the phcnomcnological 
term that was thought to count for the contribution of 
the displacement current. Comparing Eqs. ( 11 ) and ( 12 ), 



we find that the two are exactly the same when 



r,uj 
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where Ar^ 



J+^dETr(GSrgGg)(/-/) 
/ + - dETr(GSrGg)(/-/) 



(13) 

namely, the device 



is a perfect conductor with a uniform potential at 



(14) 



In other words, Eq. ( 12 ) is recovered only if the device is 



a perfect conductor with its uniform potential expressed 
by Eq. (14 1. As Eq. (11) is for the particle current, the 



fourth term in Eq. ( 12 ) is thus not from the displace- 



ment current, and is rather the correction to the particle 
current due to the induced potential of the device. It is 
important to emphasize this. 



Employing Eq. (11), we calculate the dynamic admit- 



tance of a (5, 5) carbon nanotube coupled with two alu- 
minum electrodes as depicted in Fig. la. The result is 
compared to that of TDDFT-NEGF calculation. [Uj [12] 
TDDFT-NEGF has been developed to calculate the 
time-dependent current, and employed to simulate the 
transient current through a variety of molecular and 
nanoscopic devices such as carbon nanotube based two- 
terminal device. T8J The local density approximation 
(LDA) and WBL approximation are adopted. The time- 
dependent current is evaluated, and a Fourier transform 
is performed to determine the dynamic admittance. The 
real and imaginary parts of Gll{^) are plotted in Fig. lb. 
Clearly, there is an excellent agreement between the 
frequency- and time-domain results, which confirms the 
validity of our gauge-invariant and current-continuous ac 
quantum transport theory. 

In comparison, we calculate the corresponding dy- 
namic admittances of both devices employing the phe- 
nomenological formula given in Reference [S], and the re- 
sulting GLL(a;) is plotted in Fig. lb, as well. When the 
frequency a; is zero, the phenomenological method re- 
covers both the results of the microscopic theory and 
TDDFT-NEGF. However, beyond the steady state, the 



results deviate from those of Eq. (11) and the TDDFT- 
NEGF. [11] [12] In particular, as in Fig. lb, the imaginary 
part of dynamic admittance starts to deviate immedi- 
ately at w > 0. At both low and high frequencies, 
the potential drops mostly across the CNT, which dif- 
fers drastically from the constant potential requirement 
of Eqs. (13) and (14). In fact, if the device and elec- 
trodes are mirror-symmetric, then Al = An = ^, and 
V = \Vl (with Vfj = 0). Certain constant potential dis- 
tribution Ui(r, w) = 1 would have led to the agreement 
of Eqs. (11) and (12). However, this is not the case. 



We explore how to include the contribution of the dis- 
placement current. Following the treatment of Biittiker 
and coworkers, 7^ the admittance for the displacement 
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FIG. 1. (a) The structure of a (5, 5) CNT embedded with alu- 
minum electrodes. The coordinate indicates the positions of 
atoms in real space, (b) Results from the three different meth- 
ods for an inhomogeneous Al-CNT-Al system. Solid (dashed) 
line: real (imaginary) part of the dynamic admittance from 
TDDFT-NEGF calculation; triangles (squares): real (imagi- 
nary) part from Eq. (11 1; asterisks (crosses): real (imaginary) 
part from the phenomenological formula in Reference [9]. Two 
smaller figures are the real (Re) and imaginary (Im) parts of 
both ac transport theories below 500 GHz. 
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FIG. 2. Potential distribution [uR{f,uj)] in the device region 
of the Al-CNT-Al system. 



current can be defined as 
G^M = ^ = ^^-^ = -EG«M-), (15) 



where = 0, k =^ /3. E/3 ^'^ = - E/3 Ea G^p = 



sures that the gauge invariance is satisfied. We may de- 
fine the total admittance as G^°J = Gaji -I- AqG^. Gauge 
invariance and current conservation for G*"* are satisfied 
as long as Eq = 1- This means that there are infinite 
ways to partition G^ if one is merely to ensure the gauge 
invariance and current conservation for the dynamic ad- 
mittance of the total current. Although the partition of 
G"* in Reference |9j was uniquely determined, it was based 
on the expression that violates gauge invariance for the 
dynamic admittance of the particle current. If the correct 
expression, Eq. (|9]), is used, the correction to the dynamic 
admittance would simply be zero! According to classical 
electrodynamics, the displacement current is introduced 
to ensure the current conservation, and is physically de- 
fined. According to Ampere's law, the displacement cur- 
rent should be /''(w) = ia;$(aj) , where <i>(cj) is electric 
fiux. The total current is thus J^°* — la + i^jJ^a, where 
$Q is the flux of the interface between Lead a and de- 
vice, i.e. Aa = We have thus an expression for 
the admittance from the displacement current. However, 
the calculation of the contribution from the displacement 
current can be numerically difficult, because normally the 
portions of the electrodes are included in the simulation 
box. As a result, the electric fields at the boundaries of 
the electrodes are very small, and and <i> may thus be 
very close to zero. The evaluation of Aa may encounter 
the 0^/0^ problem, which renders the accurate evalua- 
tion of Aa difficult. Moreover, Aa depends sensitively 
on the exact location of the boundary. Alternatively, 
inspired by Biittiker et al in Reference |7j, we introduce 
Terminal 0, with its dynamic admittance component ex- 
pressed as Gop{uj) — G^(cj) — —J^a^afsii^)- As a re- 
sult, the total curent conserves, since Eq=o — 

The microscopic ac quantum transport theory devel- 
oped here is gauge-invariant but not current-conserved. 
Q{uj) depends on the size of the device region. As more 
of the leads are included in the simulation box, Q{uj) be- 
comes less. In fact, if the simulation box is infinitely long 
and includes infinite portion of the leads, Q{uj) is zero; 
thus, current conservation holds. However, since the sim- 
ulation box is finite in practice, Q{u}) is usually nonzero, 
and current conservation is thus not guaranteed. 



ACKNOWLEDGMENTS 

We would like to thank Prof. Jian Wang and 
Prof. Hong Guo for helpful discussion. The fi- 
nancial support from the Hong Kong Research Grant 
Council (HKU 700808P, HKU 700909P, HKU 70071 IP, 
HKUST9/CRF/08, AoE/P-04/08), the Natural Science 
Foundation of China (Nos. 21103157 and 21233007), and 
the Fundamental Research Funds for Central Universities 
(Nos. 2340000025 and 2340000034) is gratefully acknowl- 
edged. 



5 



[1] Y. Q. Xue , S. Datta, and M. A. Ratner, |Chem. PhysT 



281, 151 (2002) 



[2] S. Datta, Nanotechnology 15, S433 (2004) 
[3] J.Tay lor, H. Guo, and J. Wang, Phys. Rev . B 63, 245407 
[ (2001) 



[4J N. D. Lang, ,Phys. Rev. B 52, 5335 (1995) 



[5] N. D. Lang and P. Avouris, [Phys. Rev. Lett. 84, 358 



(2000) 



[6j S. H. ke, H. U. Barange r, and W. T. Yang, ^J. Am. 
Chem. Soc. 126, 15897 (2004) 



[7j M. Biitti ker, A. Pretre, and H. Thomas, [Phys. Rev. Lett 



70, 4114 (1993) 



[8j M . P. Anantram and S. Datta, [Phys. Rev. B 51, 7632 



(1995) 



"WB . G. Wang, J. Wang, and H. Guo, |Phys. Rev. Lett. 82] 



398 (1999) 

[lOJ T. Yamamoto, K. Sasaoka, S. Watanabe, and K. Watan- 
abe, .Phys. Rev. B 81, 115448 (2010), 



[11] X. Zheng, F. Wang, C. Y. Yam, Y . Mo, and G. H. Chen, 

|Phys. Rev. B 75, 195127 (2007) | 
[12] X. Zheng; G. H. Chen, Y. Mo, S. K. Koo, H. Tian, C. Y. 

Yam, and Y. J. Yan, J. Che m. Phys. 133, 114101 (2010) 
[13] J. N. Zhuang, L. Zhang, and J. Wang, AIP Advances y 
P 04 2180 (2011) 

[14] Y. D. Wei and J. Wang, [Piys! Rev. B 79, 195315 (2009)] 
[15] A. P. J auho, N. S. Wingreen, and Y. Meir, |Phys. Rev. 



B 5 0, 5528 (1994) 
[16] J. Maciejko, J. Wang, and H. Guo, |Phys. Rev. B 74,| 



085324 (2006), 



[17] Y. Mo, X. Zheng, G. H . Chen, and Y. J. Yan, |J. PhysTT 



Condens. Matter 21, 355301 (2009) 
[18] C. Y. Yam, Y. Mo, F. Wang, X. B. Li, G. H. Chen, 
X. Zheng, Y. Matsuda, J. Tahir-Kheh, and W. A. God- 
dard, Nanotechnology 19, 495203 (2008), 

* Email: ghc@everest.hku.hk 



